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Metallic hydrogen and iDetalllc deuterium are both predicted to be liquids In their 
ground states over a wide range of densities. They may be described as quantisi 
liquid metals, but sinct the proton and deuteron fluids obey Ferml-Dirac and Bose- 
Elnstein statistics, respectively, their physical properties can be markedly differ- 
ent. In both there Is the likelihood of superconductivity thorugh the usual coupling 
via the density fluctuations of the Ions. But In addition the deuteronic fluid admits 
of the possibility of Bose condensation and may develop superfluid order. Finally, 
the protons may themselves pair at extremely low temperatures. 


1 . INTRODUCTION 

Under conditions of modest compression and tem- 
perature hydrogen and deuterium form diatomic- 
ally ordered Insulating states. This familiar 
situation changes radically at very high densi- 
ties where the dense phases of these systems 
undergo an Insulator-metal transition to con- 
ducting states. Though this Is expected on 
very general grounds, what cannot yet be fore- 
told with any degree of certainty Is the nature 
of the resulting high density phases, where the 
many possibilities that might be expected dif- 
fer in energy by exceedingly small amounts. 
Perhaps the most Interesting aspect of both 
metallic hydrogen and metallic deuterium Is 
that they seem to fall securely within the 
class of quantum solids and liquids [1]. Com- 
bined with their specifically metallic charac- 
ter, this places them In the somewhat rare 
category of quantum metallic solids (or as we 
shall see later, liquids). ^Thls in Itself Is 
noteworthy: however, at low temperatures 
there can be long-range ordering in both the 
electronic and ionic systems, and these too 
are Imbued with unusual characteristics. 

The purpose of this article is to discuss the 
physical attributes of some of the more physi- 
cally distinct ordered states of metallic 
hydrogen and metallic deuterium at T ■ 0 and 
nearby. In many respects the experimental 
fabrication of these metals still represents 
one of the major goals of high pressure 
physics. As noted by Ginzburg [2], the produc- 
tion of and understanding of new materials and 
substances Is a fundamental challenge. In this 
context the metallization of hydrogen has been 
and continues to be as difficult as It is fund- 
amental. The conditions required for its 
attainment involve serious problems associated 
with the ultimate strengths of materials. Thus 
the spirit of this report is very much directed 
to 'what might be' and In cataloging the remark- 
able properties below one Is nevertheless very 
conscious of the fact that the ease of theoret- 


ical conjecture bears here, apparently, an 
inverse' relationship to difficulty of experi- 
mental reality. With this cautionary remark In 
mind, we shall see tha* the inferences that can 
be made from relatively well accepted principles 
lead to some interesting predictions. These 
depend In a striking way on the quantum statis- 
tics of the systems under consideration and for 
this reason, among others no doubt, both common 
isotopes of hydrogen are worthy of detailed 
experimental examination In the high pressure 
conte.xt. 

2. THE FUNDAMENTAL HAMILTONIAN 

We begin by considering a neutral assembly of 
electrons (e) and protons or deuterons (1 * p 
or d), under conditions of density that span, 
on the one hand, the normal condensed insulating 
phases of Identifiable molecules and on the 
other high density phases exhibiting metallic 
properties. This range of densities can be con- 
veniently expressed in terms of the familiar 
linearjneasure r^, defined by Tj* O/Anpgag)^' , 
where Og is the singlo-particle density for the 
electrons averaged over the system. From the 
known density of condensed hydrogen and deuter- 
ium at 1 atmosphere, we find r • 3.07. The 
states we shall be describin', below, however, 
are characterized by r^ % 1 j. (Observe the 
experimental requirement of volume compressions 
approaching an order o' magnitude: the conven- 

ience of the linear measure of density is 
entirely theoretical.) 

Let the electi^ons in the system have a set of 
coordinates {rg}, and the protons or deutejjons 
(which we refer to as ions ) have the set {R^}. 
Then all of the phases that we will consider 
are described by a fundamental Hamiltonian 

H - H({rg).{?^}) . 0) 

All of the Interactions in the system are 
strictly Coulombic. We consider only charge 
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neutral ensembles of fixed number so that the 
thermodynamic functions are well defined. If 
Vj(k) » 4wevk^ then for a macroscopic ensemble 
of volume takes the form 

” * ^e * S ^ v^{kKN''3g(ir);g{-k)-l} (2a) 

icf^ 

* ^1 * Ifi I ''^(k){N'''p^(k)3^.(-i?)-l} (2b) 

ho 

ho 

The first two terms (2a) of H represent a well- 
known problem: they constitute the standard 
Hamiltonian of N Interactlng^electrons of mass 
nig and total kinetic energy Tg whose single 
particle density operators have Fourier compo- 
nents Pg(i(). They move In a uniform compensat- 
ing bacKg^ound at average charge density 
(eN/n). Similarly (2b) represents the Hamil- 
tonian of N protons or dguterons of mass m^, 
whose kinetic energy Is T^ and whose density 
components are pJx). They also move In a uni- 
form compensating background with average 
charge density -eN/n. With appropriate boundary 
conditions, (2a) and (2b) are themselves Hamil- 
tonians for distinct and well-defined problems. 
However, In the system under consideration they 
are coupled, and coupled strongly, by the 
attractive interaction described by (2c). In 
the thermodynamic limit, which we assime, the 
average of (2c) in states of fixed N is removed 
by the requirements of charge neutrality. 

Unlike the Hamiltonians that one might Imagine 
writing down for other simple metallic systems, 
(2) Is exact. There are no residual uncertain- 
ties stemming from the necessity of construct- 
ing a pseudopotentlal , for example. 

The Hamiltonian summarized in (2) .las consider- 
able symmetry. 

3. NORMAL GROUND STATES 

We are especially interested in ordered states 
of the system described by (2). But first we 
briefly describe the expected normal states and 
their origin. Let be the ground state wave- 
function for the electron-proton or electron- 
deuteron assembly. Then the determination of 
the ground state energy requires us to evaluate 

<f°{{?g}.(R^))|H({rg),l,v'^})|f°({?g),{R,))>. (3) 

Since (m /m )* « 0.153, and (ffig/m^)* ■ 0.128, it 
is a reasonable starting point for such a goal 
to factor the dependence on electronic and ionic 
degrees of freedom the adiabatic form; 

f°((rj.(ft;}) ♦ 'C (4) 

® ’ (R,) ® ’ 


where ij<*° Is the ground state wavjf unction for 
the Instantaneous configuration {R,} of the 
Ionic system (whose ground state, in turn. Is 
described by v’°). If the electron system 
remains in its ground state or close to It then 
It follows that the motion of the protons, or 
deuterons. Is described by the Hamiltonian 

H^({R.})-<f®° ({r,})|H({f.}.{R,})|»“ ({r J)>. 

’ {R^} ® * 

(5) 

The electronic degrees of freedom are clearly 
integrated out, but the result of tnis procedure, 
and hence of the resulting ionic motion, very 
much depends on the nature of these electronic 
states (whether Itinerant or localized, for 
example). The characteristics of these states 
deoend on density, as noted. Near r, ■ 3.1 
hydrogen and deuterium form well-defined 
diatomic molecular distance 2d(rj). (In prin- 
ciple, these units should be distinguished 
according to the states of the Ionic spins. The 
matter of ortho-para differences will not con- 
cern us further in the low density phases. ) A 
physically acceptable choice for f*” for such 
states Is a product of (N/2) two-electron func- 
tions each localized about the mean molecular 
coordinate. The electron functions In this 
choice are thus Identified by site and this 
form for therefore conflicts with the basic 
requirements of overall antisymmetry. This 
difficulty is overcome In practice by Introduc- 
ing into (5) phenomenological short range poten- 
tials which are generally taken as pairwise 
functions of the relative intermolecular separa- 
tions. At longer range, a product state based 
on an assumption of little overlap leads, as Is 
well known, to fluctuating dipole forces In 
lowest order. Thus when combined with the 
Pauli repulsion (short range) terms we arrive 
at the familiar picture of a system of interact- 
ing molecules with pairwise interactions (not 
necessarily spherically symmetric) between them. 
The ionic Hamiltonian can therefore be recast In 
the form 



nyin 


where m and m' denote molecular coordinates, and 
4 is a pair-potential . 

If the density is such that the two-particle 
functions seriously overlap, then terms beyono 
simple pair interactions displayed In (6) will 
become Important, beginning with the dipole- 
dipole-dipole term. At still higher densities 
the overlap becomes so large that the very 
starting point of the description, the assump- 
tion of well-defined localized two-electron 
states, ceases to have any physically justifia- 
ble validity. Under such conditions the density 
of electrons near the boundaries of the Wlgner- 
Seltz cells has become an appreciable fraction 
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of the average density pg, and In these circum- 
stances it is more appropriate to characterize 
the electron band states by band structure. In 
a very general s»>nse the details of the electron 
band structure must still reflect the geometri- 
cal arrangement of the protons or deuterons. 

Thus, for example, we may think of a set of 
bands which continue to display at high density 
the basic diatomic ordering so characteristic of 
the low density system. For crystalline arrange- 
ments an integral number of bands are normally 
filled, but this paradi^n of insulating config- 
urations can itself revert to metallic behavior 
through the mechanism of band overlap. That 
such overlap ran occur with diatomic ordering 
being preserved has been known for some time 
[3,4j. The correspo.iding conducting state is a 
forerunner of the more complete dissociative 
transition to a metal which, in crystals, for 
example, is expected to be represented by rela- 
tively simple structures (Bra*a1s lattices are 
the simplest). 

If metallization by band overlap is the correct 
mechanism, then the most interesting question 
concerns the density (and hence pressure) at 
which this occurs. To address this question, 
one requires a sequence of band structure cal- 
culations as a function of density for various 
separations 2d of the ions in the basis. Corre- 
spondingly one also requires a sequence of total 
energy calculations for both insulating and 
metallic phases. The metallic phases have been 
treated by a variety of methods including Wigner- 
Seitz techniques [5-9], localized orbital tech- 
niques [10-15], expansion methods [16-22] (in 
which (2c) (above) is treated as a perturbation) 
and density functional techniques [23]. The 
letter have the particular advantage that, in 
principle, they can treat metallic and molecular 
phases on the same footing with the errors of 
the method entering largely on a systematic 
basis. The combination of the results of this 


method with those of band structure calculations 
[4] leads to a prediction for band overlap metal- 
lization at around 1.8 ftoar; complete dissocia- 
tion in a static model does not appear to occur 
until much higher densities (r^ % 1.1), the 
corresponding pressure being quite dependent on 
the energies bound up in the ionic degrees of 
freedom and determined by (6). The motion of 
the ions is thus important. It is also extreme- 
ly important in determining the possible states 
of long-range order of the metallic state once 
it is formed. This motion is again determined 
by an effective io^ic Klamiltonian obtained by 
integrating out the electron degrees of freedom. 
For itinerant states appropriate to the metal 
the result is 

“ ^1 ^ ^o^^^ ^7 I ♦('^i-^i') (7) 

1,i' 

where Eq^’^s^ depending on volimie aris- 

ing from paramagnetic electron-gas energies, and 
i is a screened pair-potential for the ions. We 
note that provided r^ is held fixed, the problem 
of solving for the motion of the ions is not 
essentially different from the equivalent pro- 
blem in the insulating context. 

4. ORDERED STATES 

The description just given will determine the 
major contributions to the total energy of the 
system. The states we now discuss will 
involve ordering energies that are trivial in 
comparison. The states to be described are 
summarized in Table I. 

If hydrogen and deuterium are both in crystal- 
line states in their high density metallic 
manifestations then, with one possible excep- 
tion, it is unlikely that the fact that one is 
a Fermion system while the other a Boson will 
lead to any markrd physical distinctions. The 


System 

Normal States 

1 

Ordered States 

Crystal 

Liquid 

Crystal 

Liquid 

Electrons 

Fermi - 
Liquid 

Fermi- 

Liquid 

Super 

conductor 

Super- 

conductor 

Protons 

Quantum 

-Crystal 

Fermi - 
Liquid 

(*) 

Anisotropic- 

Superfluid 

Deuterons 

Quantum 

-Crystal 

Bose- 

Liquid 

(*) 

Bose-Condensed 

Superfluid 


Table I : States of order of electron, proton, and deuteron 

systems resulting from equation (6) and broken symmetry. The entries 
identified by' (*) might well include magnetically ordered states. 



4 


exception concerns ground state vacancies [24] 
and defects that may occur In highly quantum 
systems. The spectrum of vacancy waves might 
then be quite different. However, the presump- 
tion that the metals are Indeed crystalline, 
even In their ground states. Is Itself not at 
all well founded [1]. One can see that this 
might be the case by noting that a typical zone 
bounds. 7 phonon for the Ions will have a fre- 
quency typical of an Ionic plasma frequency, or 
an energy per Ion of about 

fiw„ % 2/?r Ry. (8) 

s 1 

For protons, for example, this amounts to 0.014 
Ry at rj • 2.0 and 0.04 Ry at r- « 1.0. These 
energies are far In excess of the energy dif- 
ferences per Ion (typically mllllflydbergs) char- 
acteristic of different static structures [22]. 
Put another way, there Is sufficient zero point 
energy In the Ionic degrees of freedom to cause 
continuous rearrangement between a variety of 
connon crystal structures. Ihough this arguixnt 
is qualitatively correct, in practice it Is nec- 
essary to treat electron and Ionic degrees of 
freedom self-consistently [25]. Nevertheless, by 
direct simulation methods L1]i It can be esta- 
blished that tne conclusion to which one Is 
being led, namely that llquid-llke ground states 
may be energetically preferred. Is in fact con- 
firmed, at least over a specified range of den- . 
si ties. If the ground state Is a liquid, how- 
ever, thr.ii the quantum statistical differences 
between hydrogen and deuterium become extremely 
important. 

4.1. Liquid Metallic Hydrogen 

Even If the density Is such that metallic hydro- 
gen Is a (quantum) crystal, the melting point Is 
expected to be very low [1]. It has been point- 
ed out that crystalline forms of metallic hydro- 
gen are likely to be superconductors with high 
transition temperatures [26-28]. Consequently, 
for crystalline states we can also Imagine that 
metallic hydrogen will be a superconductor at 
least up to Its melting point. If we now assume 
that the lack of crystallinity Is not Inimical 
to superconductive pairing (and the existence of 
amorphous or glassy superconductors certainly 
gives support to this view) then we can conclude 
that metallic hydrogen will order, even as a 
liquid. Into a superconducting state. The ques- 
tion is whether the transition temperature Is 
lower or higher than In the corresponding crystal 
and whether the standard theory can be applied. 
The answer to the latter Is that It can: to 
obtain T^ we need to solve the Ellashberg equa- 
tions for the gap function and to find the tem- 
perature at which the gap Is just suppressed by 
a pair breaking term. To carry.out this proce- 
dure, the Ellashberg function must be 
obtained from a description of density fluctua- 
tions which Is more general than that used for 
crystals. In turn, this requires some knowledge 
of the frequency dependent density-density 


response function x(Ri<u) whose Imaginary part 
appears explicitly [29] in the Ellashberg func- 
tion. This can be obtained from memory func- 
tion techniques, for example, but Is most easily 
given by a generalization of the RPA like form: 




ILfjSUii] 

|-T(<l.u)x°(<1fi*i) 


^9) 


where response function for non- 

interacting protons and f(q,u) is the so-called 
polarization potential [30]. It can be obtained 
from the reasonable assertion that It is a con- 
tinuous function with the appearance of a 
screened Coulomb Interaction at long range, and 
taking a constant value at short range. The 
constant Is fixed by the calculated sound speed 
as q ♦0. From (9) and from the electron-proton 
Interaction v^Wq) we obtain the Ellashberg 
function according to 

a^F(u) « N(0) I ^dq(q/2kp)|Vg^(q)^x‘'(q.u>)/’f- 
■^0 

( 10 ) 


This Is the required Input for the gap function 
b(;j): the solution of the Ellashberg equations 
[31] Is relatively straightforward. The result 
Ing transition temperature Is 'u 10‘°K and Is a 
fairly strong function of density [29]. These 
temperatures do Indeed exceed the melting tem- 
perature in certain density ranges so that the 
assumption that a liquid superconducting state 
exists appears to be substantiated. 


It Is worth noting that If the density is such 
that metallic hydrogen Is in a liquid state, 
but the temperature Is such that the system Is 
normal, the properties of this 'normal' state 
are still rather remarkable. For It then con- 
stitutes a two-component Fermi liquid with one 
component (the electrons) In the high density 
limit and the other (the protons) In the low 
density limit. Both Fermi liquids have Identi- 
cal Fermi surfaces but very different Fermi 
temperatures (Tp„ -v 10*3 Tp^). From the stand- 
point of experimental detection of a(statlc) 
pressure generated phase of metallic hydrogen. 
It Is also Important to understand the proper- 
ties of this phase. This can be achieved most 
easily through a straightforward generalization 
of Landau Fermi liquid theory to two (charged) 
components, both for equilibrium properties 
[32] and transport properties [’3]. One of the 
most striking equilibrium properties is the 


appearance of a giant linear specific heat at 
low temperatures? ft Is a manifestation of the 



phonon behavior expected when the system freezes. 

The transport properties are no less dramatic. 

To study these one begins with a two-component 
general Ization [33] of the Landau Silin trans- 
port equations describing the momentum, space, 
and time development of the electron and 
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proton quasi -particle distribution functions. 

It is possible to obtain closed form expressions 
for the electrical conductivity, thermal conduc- 
tivity <. viscosity n. *nd spin diffusion coef- 
ficients. As T ♦ 0 all of these quantities 
diverge because of the phase-space considerations 
controlling the scattering. The electrical 
resistivity, in particular, vanishes as T‘, 
again in marked contrast to the solid where, in 
the absence of impurities, it would vanish as 
T*. Most importantly, the transport coeffi- 
cients of the liquid are found to be at least a 
factor {mg/ra_)2 smaller than In the correspond- 
ing solid metal (and in the case of the resis- 
tivity, the inverse of the conductivity, the 
factor can approach 10®!). 

Finally, in a fluid phase of metallic hydrogen, 
the Fermi temperature of the protons <s 
(m./m_)Tr, where Tc„ is the Fermi temperaf v' 
of t^e felectrons. For rj 'v- 1-2, Tp- 10‘®K. 

As noted above, there are effective' (screened) 
interactions $ bet.«en the protons which have 
a repulsive region at short range, and an 
attractive and oscillatory region at long range 
[1]. Given this form of the interaction and 
with the example of He^ before us, it is possi- 
ble that at extremely low temperatures the 
protons may themselves ue unstable toward pair- 
ing. The temperature at which such a protonic 
superfluid will form will be in the neighborhood 
of 10'3 Tpp. Whether the potential is suffi- 
ciently repulsive at short range to favor p-wave 
over s-wave pairing has yet to be determined. 

4.2. Liquid Metallic Deuterium 

A possible low temperature liquid metallic phase 
of deuterium represents a system with “mixed 
statistics," the electrons constituting a Fermi 
liquid, as before, but the deuterons belonging 
to the class of spin-1 Bosons. Unlike helium 
(spin-zero Bosons) the presence of the addition- 
al spin degree of freedom leads to a new branch 
in the quasi-particle spectrum even though 
there are no explicitly spin-dependent terms in 
the Hamiltonian (2). The dispersion of the branch 
follows a k^ behavior at small k [35], and 
therefore differs in an essential way from the 
linear type of dispersions so characteristic of 
phonon-like quasi-particle excitations. There 
are thus 3 separate types of elementary excita- 
tions: the "phonons," the spin excitations, and 
electron-hole pairs. The system is quite 
unusual [36] and certainly contrasts markedly 
with, say, He’-He’ mixtures. 

At low temperatures we may certainly expect, as 
with metallic hydrogen, that electron pairing 
can occur through the density fluctuations of 
th{ ions. The crucial difference is that 
x(q,ui) (see equation (9)) now reflects the Boson 
character of the ionic system. Again the 
approach is very similar: x is constructed with 
the polarization-potential model [30] and the 
Eliashberg equations again solved. The resulting 
transition temperatures are somewhat lower than 


in metallic hydrogen [37], but still in the 
neighborhood of 10‘K. 

At even lower temperatures a potentially more 
interesting ordering can occur in liquid metal- 
lic deuterium. This is a superfluid transition 
associated with the Bose-condensation anticipat- 
ed for the spin-1 system. At rj ■ 1.5 an 
estimate based on non-interacting Bosons gives 
Tgr Hi 30“K for this system. This suggests the 
rather intriguing possibility that we have a 
charged system in which the electrons display 
superconducting order while at the same time 
the deuterons display superfluid order. If it 
exists it is difficult to imagine a more remark- 
able state of condensed matter. At the very 
least it raises a mmber of Interesting theoret- 
ical issues, for example, whether or not it is 
valid to describe the system with a two-fluid 
hydrodynamical picture, and whether or not per- 
sistent deuteron currents will actually exist, 
that is to say, whether combined superfluidity 
and superconductivity is stable. 

The transport properties of the system are 
expected to be unusually r-*wn in their possible 
temperature dependences. This is because in 
setting up the appropriate kinetic equations for 
the three types of quasi -particles it is recog- 
nized that their associated statistics have con- 
siderable bearing on the standard phase-space 
arguments determining the scattering character- 
istics. To give but one example, if the tem- 
perature exceeds both Tdq for the deuterons and 
T(; for the electrons, then the low temperature 
behavior of the electrical resistivity is expect- 
ed to be linear in the temperature [36] rather 
than quadratic as is the case for the normal 
electron-proton fluid. The reason is that the 
phase space limitation on the scattering cross- 
section Imposed by the Fermi-Dirac distribution 
of the protons is essentially removed when the 
protons are replaced by deuterons. A further 
consequence is that the magnitude of the low 
temperature resistivity of normal liquid metal- 
lic deuterium should greatly exceed that of 
liquid metallic hydrogen. As T falls below Tg/. 
the deuterons will Bose condense and this tran- 
sition should be manifested in the electrical 
transport properties, both in critical behavior 
near Tg^ itself and in the form of the transport 
coefficients well below Tgj- [36]. In particular, 
the manner ii' which electrons cause excitations 
between the condensate (presumably in one sub- 
space of the spins) and the excitation gas (in 
the remaining subspaces) should be particularly 
interesting. 

It was remarked above that the Hamiltonian (2) 
has considerable symmetry. In the macroscopic 
manifestations of (2) this symmetry is broken, 
the resulting states possessing physical proper- 
ties that may be unique among the elements. The 
states we have discussed are probably onl> the 
simplest that can be imagined. There may be 
more complex possibilities involving magnetic 
order. 
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